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The optical effects of the light produced by lightning axe of interest to atmospheric 
scientists for a number of reasons (4,7,13). We mention two techniques used to explain the 
nature of these effects: Monte Carlo simulation (13) and an equivalent medium approach 
(10). In the Monte Carlo approach, paths of individual photons axe simulated; a photon 
is said to be scattered if it escapes the cloud, otherwise it is absorbed. In the equivalent 
medium approach, the cloud is replaced by a single obstacle whose properties axe specified 
by bulk parameters obtained by methods due to Twersky. See Phanord (10) for references. 

In this report, we use Boltzmann transport theory to obtain photon intensities. The 
photons are treated like a Lorentz gas. We consider only elastic scattering and neglect 
gravitational effects. Water droplets comprising a cuboidal cloud are assumed to be 
spherical and homogeneous. Furthermore, we assume that the distribution of droplets in 
the cloud is uniform and that scattering by air molecules is negligible. The time dependence 
and five dimensional nature of this problem make it particularly difficult; neither analytic 
nor numerical solutions axe known (3). 

We begin with the single speed Boltzmann transport equation (3) 


^ = -cCl ■ VI 4- —Ku o f Vldtf - Kcl + cs, I = I (r, Q, t ) , 
dt 47 t J v ' 


(n), 
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and 


where I is the photon intensity which depends on position (r), velocity direction 
time ( t ). Inside the integral, I is considered a function of Cl 1 . Here, c is the speed of light, 
K = i/A, where A is the mean free path, u 0 is the single scattering albedo, V= v(Cl’ • Cl) 
is the scattering phase function, and s is a source term. The intensity in a volume V 
bounded by a surface S is uniquely determined by the initial intensity in V , the sources 
in V, and the intensity incident on S (1). 

In order to reduce the number of dimensions, we use the Pn approximation (1,3). The 
intensity, phase function, and source term are all expanded in series of spherical harmonics 

= pficosdy™*-, Ci = Ci{9 , <t>y, 


[2] 


/(r, n, t) = £ 6 :c(f>). k = c(r, ‘).E=EE' 

n=0 m=— n 
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The addition theorem for spherical harmonics (12) allows us to write 

p(sV ■ !i) = £ ( - (sV) (n) . 

In particular, if we use the Henyey- Greenstein function (14), d n = (2n + l)<7 n j where g 
is the asymmetry factor. The integral in [1] may be evaluated using orthogonality. The 
advection term of contains an inner product of Cl and a gradient which we may write as 


[3] 
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We absorb the trigonometric functions into the spherical harmonics comprising I with the 
aid of the recurrsion relations for the Legendre functions (9) and orthogonality is used to 
obtain the coupled set of partial differential equations, 
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c dt 
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where any function b™ with m > n or n < 0 is identically zero. 

We work with a truncated version of [5] motivated by the Eddington approximation 
(2), b™= 0 for 77 > 2 . Using the notation b 0 o = / 1? b\ — = / 2 , l — A’ = A’ 

we have 


id/, 
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s = s„ 


where we have assumed that s™=o, 77 > 1 . This truncation does not constitute a consistent 
solution to the transport equation, a fact implicit in (2,5,6) and explicit in (1). Orthogo- 
nality may be used to produce an additional five equations which are neglected. 

The conditions at an interface between a convex cloud and a vacuum are given by (1,3) 


l(p, U, = 0, (l • 77 < 0, [7] 

where p is a position vector on the cloud’s boundary and n is the outward normal. 
Physically, this condition corresponds to the requirement that for convex cloud geometries, 
photons which escape the cloud cannot reenter. With the simplified representation given 
by the Eddington approximation, it is generally impossible to satisfy [7] exactly. The 
condition usually used is the Marshak boundary condition (1,2, 3, 5, 6), 

F n = - j Sl-hl(p, Cl, f)dft = 0 , [ 8 ] 

(th< 0 


where F n is the irradiance in the direction h. The integral of F n over a portion of the 
interface yields the number of photons reentering the cloud through that portion. 
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Consider a plane surface separating the cloud and the region outside with normal h. 
Since we may take the area to be vanishingly small and still require that no photons reenter 
through this surface, [8] may be rewritten as 


2X TT 

/<#„/ d0 n sin 0 n cos 6 n I[9,<f>]d£l = 0, [9] 

O * 

a 

where I[0 , d>] is evaluated on the cloud boundary and 6 n is measured from h considered as 
a polar axis. An appropriate rotation of coordinates yields the condition 

fi = -{n x f a +n y f 3 +n z f A ), h = n x i+n x y+n x z [10] 

3 

on an interface. This form is valid for any convex shape and is a generalization of a similar 
condition given by Davies (3). 

We use Laplace transforms in [6] and obtain 

& 2 Fl = V 3 F X + T, k* = i(p o K+-^(pj<+-^), Ts 3 (A^)(^7 +S )- 

where s is the transform variable, Fi and S are the Laplace transforms of /i and s 
respectively, /i(o) is the initial value of f lt and we have assumed that / 2 = fs = f\ — 0 
at t = 0. The boundary condition becomes 


f^F 1= 0, h = ^(-+PJ<). [12] 

on 2 \ c / 

We consider a cuboidal cloud centered at the origin with boundaries coinciding with 
x _ y _ ±1^ z — ±1, Equations [11,12] are solved using finite Fourier transforms 

(14) to obtain 


f.-E 
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where T is the finite Fourier transform of T and £p, f]q , v r ^re the roots of 


c* h P h 7 h 

tan 2^ = v tan 2 ,i = v Un r r= v, 


[14] 


with positive real parts. 

Due to the transcendental nature of the functions defined by [14], it is necessary to 
resort to numerical methods. Fortunately, it is possible to approximately invert Laplace 
transforms with a knowledge of Fl (s), 5 = 1,2,3 ,... • Complete details are provided in 
Lanczos (8). 
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Once f 1 is known, / 2 , / 3 , / 4 may be found using [6]. The photon intensity at a point 
outside of the cloud in a direction tl may be traced back to its value on the cloud’s surface 
at an earlier time t' (1), i.e., 7^r, tj = I (r — c(l(t — t'),Cl, t'^, where t 1 is chosen so 

that r — c(l(t — t 1 ) = p. 

We are presently working on implementing the algorithms described above and on a 
numerical procedure which is not limited to the Eddington approximation. We intend to 
generalize the results to other cloud shapes and compare results with existing Monte Carlo 
simulations as well as with results using an equivalent medium approach. 
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